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摘 要 : 精准 施肥 是 减少 农业 面 源 污染 的 重要 技术 之 一 , 而 土壤 养分 测试 与 作物 营养 诊断 是 其 实施 的 技术 保障 ,特别 是 在 农 
业 规 模 化 经 营 方式 下 , 急需 发 展 快速 、 经 济 、 无 损 的 作物 氮 素 营养 诊断 技术 。 本 文 在 应 用 数字 图 像 进 行 冬小麦 、 夏 玉米 氮 素 
营养 诊断 研究 的 基础 上 , 将 数码 相机 搭载 到 无 人 机 上 ，, 利用 无 人 机 航拍 技术 采集 作物 冠 层 数字 图 像 ,研究 不 同 航拍 高 度 下 完 
层 图 像 相关 色彩 参数 反 演 冬 小 麦 和 夏 玉米 氮 素 营养 状态 的 差异 ， 以 确定 适宜 的 航拍 高 度 与 敏感 的 色彩 参数 ,建立 利用 无 人 机 
航拍 数字 图 像 诊断 冬小麦 和 夏 玉 米 氮 素 营养 状态 的 模型 。 研究 结果 表明 : 在 冬小麦 拔节 期 适宜 的 航拍 高 度 是 16 m, 敏感 的 色 
彩 参数 是 可 见 光大 气 阻抗 植被 指数 (VARD，, 诊断 模型 为 : 冬小麦 茎 基部 硝酸 盐 浓 度 =2.103 4el8874vARI;， 夏 玉米 大 喇叭 口 期 适 
宜 的 航拍 高 度 是 50 m， 敏 感 色彩 参数 是 蓝光 标准 化 值 [BRHGHB)],， 诊断 模型 为 : 夏 玉米 第 1 完全 展开 叶 叶 脉 硝酸 盐 浓 度 
=1.526.1032.[B/R+G+B)]5045。 依 据 建立 的 航拍 方法 与 诊断 模型 ， 分 别 对 冬小麦 、 夏 玉米 进行 了 氮 素 状态 监测 的 验证 ,结果 表 
明 诊 断 结果 与 冬小麦 、 夏 玉米 实测 数据 的 决定 系数 分 别 为 0.80 和 0.85， 目 均 在 P<0.01 水 平 显著 相关 。 最 后 将 研究 结果 进行 
应 用 ,生成 了 冬小麦 、 夏 玉米 氮肥 追肥 作业 图 。 利 用 无 人 机 搭载 数码 相机 对 冬小麦 、 夏 玉米 进行 氮 素 营养 诊断 简单 、 可 行 , 但 
仍 有 一 些 技术 细节 需要 完善 ， 以 提高 该 技术 的 实用 性 。 

关键 词 : 精准 施肥 ; 无 人 机 ; 数码 相机 ; 色彩 参数 ; 航拍 高 度 ; 硝酸 盐 浓 度 ; 氮 素 营养 诊断 ; 冬小麦 ; 夏 玉米 
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Diagnosis of nitrogen nutrition of winter wheat and summer corn using images 
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from digital camera equipped on unmanned aerial vehicle 
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Abstract: The large number of population and scarce land resources require crop yield per unit area to be increased greatly to ensure 
the food supplies and food security in China. To meet this requirement, chemical fertilizer application ls an important and simple 
measure. Consequently, the agricultural non-point pollution caused by chemical fertilizer is becoming more and more serious. 
Precision fertilization is one of the important technologies to control agricultural non-point source pollution, while soil nutrient 
testing and crop nutrition diagnosis are indispensable works of precision fertilization. With the increasing of agriculture scale 
management, it is urgent to develop fast, economic and nondestructive techniques of the nitrogen nutrition diagnosis of crops. In this 
study, in the light of related research of diagnosis of nitrogen nutrition in winter wheat and summer corn using digital photographs, 
cost-effective unmanned aerial vehicle (UAV) system was introduced and digital camera mounted on its platform was used to obtain 
the canopy photographs of crops. To determine the appropriate aerial heights and sensitive color parameters of aerial photographs for 


the diagnosis of nitrogen nutrition of winter wheat and summer corn, the relationships between a variety of color parameters retrieved 
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from the aerial photographs obtained at different heights and the status of nitrogen nutrition of winter wheat and summer corn 
measured in situ in 2015 were analyzed. It was indicated that VARI (visible atmospherically resistant index based on digital 
photography) and B/(R+G+B) were the most sensitive color parameters for the diagnosis of nitrogen nutrition in winter wheat and 
summer corn, respectively. For the winter wheat, 16 m was the suitable aerial height at jointing stage, and the diagnosis model was: 


ISSTVARL, 50 m was the best aerial height for summer corn at big trumpet stage, and the 


stem sap nitrate concentration = 2.103 4e 
diagnosis model was: Vein nitrate concentration = 1.526:10™: [B/(R+G+B)]0 5 The established diagnose models were verified at 
the same field in 2016. The results showed that the determination coefficients of relationship between predicted and observed nitrate 
nitrogen concentrations of stem sap of winter wheat and vein of summer corn were 0.80 and 0.85, respectively, with significant levels 
of P<0.01.1n 2017, the method was applied to the diagnosis of nitrogen nutrition of winter wheat and summer corn in Nanpi County, 
Hebei Province, and the prescription maps of top-dressing nitrogen fertilizer for winter wheat and summer corn were made. The 
nitrogen fertilization rates of winter wheat and summer corn given by the prescription maps were fitting well with the growth status 
of crops. As shown in this study, diagnosis of nitrogen nutrition of winter wheat and summer corn using photographs from digital 
camera equipped on UAV was a cost-effective way to achieve precision management of fertilization. Beyond the achievements of 
this study, there were still problems should be solved to enhance the practicability and applicability of this method. 

Keywords: Precision fertilization; Unmanned aerial vehicle (UAV); Digital camera; Color parameters; Aerial height; Nitrate nitrogen 


concentration; Diagnosis of nitrogen nutrition; Winter wheat; Summer corn 


我 国人 多 地 少 , 要 保证 粮食 供应 必须 提高 作物 单产 ,而 施用 化 肥 是 提高 作物 产量 的 重要 措施 。 近 年 来 
我 国 粮食 产量 的 连续 增长 , 在 很 大 程度 上 依赖 于 化 肥 的 施用 帆 。 较 高 的 化 肥 用 量 和 较 低 的 养分 利用 率 已 经 
对 环境 产生 了 不 良 影响 , 粮食 安全 政策 与 农业 面 源 污染 的 关系 得 到 了 广泛 关注 , 为 此 ,农业 部 提出 了 《到 
2020 年 化 肥 使 用 量 零增长 行动 方案 》, 其 中 精准 施肥 是 重要 技术 之 一 站。 实施 精准 施肥 的 前 提 是 土壤 养分 
测试 与 作物 营养 诊断 。 土 壤 养 分 测试 技术 相对 成 熟 ， 农业 部 在 全 国 范 围 内 组 织 开展 了 测 土 配方 施肥 普及 行 
动 趾 ,而 在 作物 氮 素 营养 诊断 方面 , 传统 的 实验 室 方法 虽然 精准 但 费时 费力 , 难以 满足 作物 快速 生长 时 期 农 
田 氮 肥 管理 的 需求 。 因 此 ,发 展 快速 、 经 济 、 无 损 的 作物 氮 素 营养 诊断 方法 对 于 实现 合理 施用 氮肥 、 提 升 
作物 品质 、 增 加 农民 收入 和 减少 面 源 污染 都 具有 重要 的 应 用 价值 。 

作物 氮 素 无 损 诊断 随 着 相关 技术 的 提高 由 定性 或 半 定 量 逐 渐 癌 精准 定量 的 方向 发 展 , 主要 包括 叶绿素 
仪 、 腺 感 技 术 、 荧 光 诊 断 法 和 冠 层 色彩 分 析 法 等 中。 其 中 , 冠 层 色 彩 分析 法 基于 作物 缺 氮 时 叶 色 的 变化 来 
诊断 其 氮 素 营养 状况 。《 沈 氏 农 书 》 中 就 有 关于 根据 水 称 (Oryza sativa) 叶 色 追 施 穗 肥 的 记载 。20 世纪 60 年 
代 初 ， 陶 勤 南 等 中 通过 比 色 卡 的 方法 评价 水 稳 叶 色 并 对 其 黑 黄 变化 的 生理 基础 做 了 系列 研究 。 近 期 随 着 数 
字 照 相 技 术 的 发 展 ， 冠 层 色彩 分 析 法 进入 精准 定量 研究 阶段 ,并 在 多 种 作物 上 得 到 应 用 。 张 立 周 等 "| 基于 
养分 定位 试验 研究 了 小 麦 (TYiticum aestivum)、 玉 米 (Zea mays) 冠 层 数 字 图 像 色彩 参数 与 其 传统 氮 素 营养 参数 
之 间 的 关系 ,筛选 出 绿 光 标准 化 值 和 蓝光 标准 化 值 分 别 用 于 小 麦 和 玉米 的 氮 素 营养 诊断 。 王 连 君 等 四 发 现 
绿 光 /蓝光 可 用 于 葡萄 (iis vinifera) 氮 素 营 养 的 诊断 。 魏 全 全 等 中 采用 冬 油菜 (Brassica campestris) 冠 层 数 字 
图 像 的 红 光 标准 化 值 对 其 氮 素 营养 进行 诊断 。 另 外 , 数字 图 像 氮 素 诊断 技术 在 马 铃 划 (So1arzm 
tpberosz1n)lg、 草 莓 (Fragaria ananassaU、 棉 花 (Gossypizm spp.) "等 多 种 作物 上 均 得 到 研究 和 应 用 ， 充 分 
显示 该 方法 可 以 作为 优化 配方 施肥 的 配套 技术 辅助 实现 氮肥 的 精准 管理 。 

随 着 农村 土地 流转 量 的 增加 ,农业 适度 规模 经 营 成 为 新 的 趋势 并 在 不 断 扩 大 , 采用 数字 图 像 氮 素 诊断 
技术 进行 氮肥 精准 管理 必须 解决 快速 获取 大 范围 作物 冠 层 图 像 的 问题 。 无 人 机 技术 的 发 展 为 低空 匠 感 提供 
了 可 行 的 搭载 平台 。 高 林 等 "利用 数码 相机 与 旋 由 无 人 机 组 成 的 壳 感 监测 系统 对 冬小麦 叶 面积 指数 进行 了 
确 估算 。 祝 锦 霞 等 ("采用 旋翼 式 无 人 航空 摄影 平台 对 水 稻 氮 素 营 养 诊 断 开 展 了 研究 。 由 于 无 人 机 遥感 的 
时 效 、 高 分 辨 率 和 低 成 本 的 特点 , 在 田 块 尺度 农 情 监测 与 决策 方面 具有 很 好 的 应 用 前 景 。 低 空 对 作物 冠 
航拍 不 同 于 地 面 拍照 , 在 航拍 高 度 、 速 度 、 数 码 相 机 拍照 参数 的 确定 上 要 力求 图 像 清 晰 、 拍 照 区 域 履 盖 
整 日 不 过 于 重 营 、 工 作 效 率 高 效 , 在 数字 图 像 获取 的 航拍 技术 、 氮 素 诊 断 指标 筛选 及 图 像 处 理 技 术 上 需 
进一步 的 研究 以 满足 实际 生产 的 需求 ， 而 这 方面 的 研究 还 很 欠缺 。 

基于 上 述 的 研究 基础 和 技术 需求 ， 本文 利用 无 人 机 搭载 数码 相机 对 冬小麦 、 夏 玉米 冠 层 进行 航拍 , 结 
合作 物 常 规 氮 素 分 析 ， 对 航拍 数字 图 像 进行 分 析 以 确定 适宜 的 航拍 高 度 ， 和 划 选 出 敏感 的 色彩 参数 ， 实 现 对 
作物 氮 素 营养 状态 的 精准 诊断 ， 以 期 为 规模 化 农业 经 营 方式 下 农田 氮肥 的 精准 管理 提供 技术 支持 。 
1 材料 与 方法 
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1.1 大 田 试验 设计 

为 了 监测 不 同 程度 氮 素 短缺 在 冬小麦 、 夏 玉米 叶 色 上 的 反应 差异 , 试验 在 长 期 肥料 定位 试验 地 中 进行 
5。 试验 地 位 于 中 国 科学 院 栾 城 农业 生态 系统 试验 站 (37.88°N, 114.68°E), 该 站 属 暖 温带 半 湿 润 季 风气 候 , 
代表 华北 平原 北部 典型 湖 褐 土 高 产 农 业 生 态 类 型 ,种 植 制度 以 冬小麦 - 夏 玉米 轮作 为 主 。 肥 料 定位 试验 自 
1997 年 开始 一 直 延 续 至 今 , 冬小麦 种 植 品 种 为 “ 科 农 199”, 夏 玉 米 品 种 为 “ 郑 单 958”,，, 共 设 4 个 氮肥 水 
平 : 0、200 kg(N):hm .a 1 、400 kg 人 -hm a 和 600 kg(N):hm a!; 3 个 磷肥 水 平 : 0、32.5 kg(P).hm 2.a 1 
和 64 kg(P).hm a 1; 2 个 钾肥 水 平 : 0 和 150 kg(K):hm “a 1。 试 验 采用 正 交 不 完全 设计 共 16 个 处 理 , 每 个 
处 理 3 次 重复 共计 48 个 样 区 。 为 了 避免 磷肥 严重 短缺 的 影响 , 本 研究 选择 施用 磷肥 的 处 理 作 为 供 试 样 地 ， 
共 33 个 样 区 参与 本 试验 的 采样 与 分 析 。 
1.2 无 人 机 搭载 相机 与 航拍 方法 

本 研究 采用 大 疆 S900 六 悬臂 无 人 机 搭载 尼康 D7000 数码 相机 进行 冬小麦 和 夏 玉 米 冠 层 的 航拍 。 研 究 
组 自制 了 减 震 平台 将 数码 相机 搭载 到 无 人 机 上 ,搭载 平台 如 图 1 所 示 。 为 了 减少 无 人 机 振动 对 成 像 效果 的 
影响 ,在 连接 杆 与 无 人 机 及 平台 的 连接 处 均 安装 了 减 震 球 ， 以 减少 了 振动 传递 。 搭 载 平台 均 采 用 碳纤维 材 
料 , 保证 结构 强度 同时 减少 重量 ,避免 额外 负担 降低 无 人 机 飞行 时 间 。 相 机 水 平 向 下 檬 入 搭载 平台 ,航拍 时 
将 其 固定 ,避免 滑 动 或 掉 落 。 由 于 相机 不 具备 GPS 模块 ， 需 要 配备 专用 GPS 接收 器 , 在 航拍 时 为 相机 提供 
经 纬度 位 置信 息 ， 并 将 其 记录 在 图 像 的 EXIF 信息 中 ,以 备 后 期 图 像 处 理 时 需要 。 
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图 1 无 人 机 搭载 数码 相机 平台 结构 图 


Fig.1 Structure diagram ofunmanned aerial vehicle platform equipped with digital camera 


无 人 机 航拍 时 利用 地 面 站 软件 规划 好 起 落 点 、 飞 行 轨迹 、 高 度 、 上 升 和 水 平 飞行 速度 等 参数 ， 然 后 上 
传 到 无 人 机 飞 控 ， 实 现在 固定 航线 和 高 度 上 自动 航拍 。 为 了 保持 相机 尽量 水 平 获 得 正 射 图 像 ， 将 水 平 飞行 
速度 设 定 为 3 ms-'。 航拍 选择 在 晴朗 无 云 、 风 力 较 小 的 天 气 , 时间 在 11:00 一 13:00, 避免 云 和 太阳 光 强 变化 
对 图 像 的 有 影响。 尼康 D7000 数码 相机 配备 了 18~300 mm 1V3.5-5.6G 镜头 , 航拍 时 将 焦距 设 定 为 18 mm， 以 
JPG 格式 保存 图 像 。 

本 研究 分 别 于 2015 年 和 2016 年 在 肥料 定位 试验 地 进行 航拍 , 2017 年 在 南 皮 县 进行 了 应 用 。 其 中 , 2015 
年 数据 用 于 筛选 色彩 参数 建立 诊断 指标 , 2016 年 用 于 指标 的 验证 。 航 拍 时 间 选 择 在 冬小麦 拔节 期 (2015 年 4 
月 9 日 和 2016 年 4 月 8 日 ) 和 夏 玉米 大 喇叭 口 期 (2015 年 7 月 28 日 和 2016 年 7 月 22 日 )。 为 了 研究 不 同 高 
度 下 航拍 图 像 对 冬小麦 、 夏 玉米 氮 素 状态 反 演 的 差异 , 冬小麦 航拍 高 度 选择 8 m、10 m、12 m、14 m、16 m 
和 55 m， 夏 玉米 航拍 高 度 选择 10 m、15 m、20 m 和 50m。 航 拍 时 相机 采用 2s 间隔 进行 自动 拍照 ,同时 在 
定位 试验 设 地 置 数字 标牌 ， 以 便 将 图 像 与 地 面 采 样 数据 进 行 对 应 。2017 年 4 月 14 日 和 7 月 27 日, 在 河北 
省 南 皮 县 东 五 拨 村 (38.068*N, 116.825"E) 分 别 对 冬小麦 和 夏 玉 米 地 进行 航拍 , 利用 本 研究 建立 的 诊断 方法 进 
行 氮肥 追肥 推荐 ， 指 导 农 户 精 准 施 肥 。 
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对 33 个 样 区 在 航拍 时 间 步 测量 作物 氨 素 含量 。 利 
中 部 ， 每 个 样 区 测量 30 个 叶片 ， 取 其 平均 值 作为 该 区 


一 


1.4 和 冠 层 图 像 处 理 与 参数 分 析 方 法 


用 SPAD-502 叶绿素 仪 测 量 冬小麦 第 1 完全 展开 叶 的 


SPAD 数值 。 采 用 反射 仪 和 硝酸 盐 试纸 测量 冬小麦 茎 
基部 硝酸 盐 浓度 (以 下 简称 冬小麦 茎 基部 酸 盐 浓度 )。 同 样 取 夏 玉米 第 1 完全 展开 叶 中 部 叶脉 进行 挤 压 汁液 ， 
采用 同样 的 方法 测量 叶脉 硝酸 盐 浓度 (以 下 简称 玉米 叶脉 硝酸 盐 浓度 )。 


冠 层 图 像 通过 红 (R)、 绿 (G)、 蓝 (B)3 个 色彩 通道 来 反映 作物 叶片 对 光 的 反射 特性 ,利用 无 人 机 搭载 相 
机 对 作物 进行 航拍 可 以 最 大 限度 地 获得 顶部 展开 叶 的 信息 。 但 是 , 由 于 光线 遮挡 等 原因 ， 图 像 中 存在 阴影 
裸 地 、 枯 叶 等 非 作 物 叶 片 信息 ， 如 果 全 幅 图 像 均 参 与 色彩 参数 分 析 , 混入 的 干扰 信息 将 影响 到 敏感 参数 的 


时 人/ 
呆 y 


选择 ,因此 需要 通过 一 定 的 过 滤 通 道 将 非 作 物 叶 片 的 像素 剔除 ， 以 提高 敏感 色彩 参数 的 筛选 及 诊断 精度 。 


由 R、G、B 色彩 通道 相互 组 合 可 以 衍生 出 多 种 色彩 参数 ， 如 果 都 参与 同 常规 氮 素 测量 参数 的 相关 性 分 


析 ，, 工作 量 将 是 很 大 。 参 考 相关 研究 ,同时 考虑 到 比值 型 色彩 参数 对 环境 光线 变化 的 消减 作用 ， 本 研 
究 通 过 初步 筛选 ,选择 红 光 标准 化 值 [R/(R+G+B)]、 绿 光标 准 化 值 [G/(R+G+B)]、 蓝 光标 准 化 值 [B/(R+G+B)] 
以 及 可 见 光 大 气 阻抗 植被 指数 [(G-R)/A(G+R-B)，VART] 作 为 备 选 指标 , 与 同步 测量 的 作物 氮 素 指标 进行 相 
关 性 统计 分 析 ， 以 确定 理想 的 航拍 高 度 与 诊断 色彩 参数 。 


2 结果 与 分 析 
2.1 非 作物 冠 层 像素 剔除 


航拍 图 像 中 非 冬 小 麦 、 夏 玉米 叶片 像素 的 存在 会 对 参数 选择 及 诊断 造成 干扰 ,分 析 时 必须 剔除 ， 以 保 
证 只 有 作物 叶片 的 色彩 参数 参与 统计 与 分 析 。 通 过 对 冬小麦 、 夏 玉米 航拍 图 像 中 叶片 及 其 他 干扰 地 物色 彩 
参数 的 分 析 , 选择 G/(R+G+B)、R/(R+G+B) 以 及 R 共 3 个 色彩 参数 通过 阔 值 设 定 对 图 像 进 行 过 滤 。 过 滤 掉 
非 叶片 像素 的 规则 是 : G/R+G+B)>0.37 AND R/(R+G+B)<0.42 AND R>20, 其 中 G/R+G+B) 与 RA(R+G 


用 于 叶片 像素 的 选择 , R>20 用 于 过 滤 阴 影 。 图 2 为 冬小麦 、 夏 玉米 航拍 图 像 与 过 渡 后 形成 的 掩 膜 图 像 对 
可 见 航 拍 图 像 中 的 裸 地 、 标 牌 、 杂 物 、 阴 影 等 都 已 经 被 剔除 ， 这 将 提高 后 期 图 像 色 彩 参数 分 析 的 准确 性 。 


图 2 冬小麦 和 夏 玉米 航拍 图 像 中 非 叶片 像素 剔除 效果 


Fig.2 Comparison between aerial photographs of winter wheat and summer corn and their mask images 


HB) 
比 ， 


利用 过 滤 非 叶片 像素 生成 的 叶片 掩 膜 文件 (0 和 1 数值 ) 和 每 个 叶片 像素 的 3 个 色彩 通道 数据 , 计算 获得 
各 像素 的 R/(R+G+B)、G/(R+G+B)、B/(R+G+B) 和 VARI。 根 据 地 面 标牌 号 码 识别 每 个 样 区 的 航拍 图 像 ， 


算 并 统计 出 每 个 样 区 所 有 叶片 的 色彩 参数 获得 平均 值 ， 


析 ， 从 中 筛选 适宜 的 航拍 高 度 和 敏感 的 色彩 参数 ， 建 立 诊断 方程 。 


计 


与 同步 测量 的 冬小麦 、 夏 玉米 氮 素 状态 信息 进行 分 


2.2 ”冬小麦 氮 素 营养 诊断 分 析 
2.2.1 敏感 参数 选择 

手持 式 叶 绿 素 仪 (SPAD-502) 能 够 快速 、 简 单 、 无 损 地 测量 植物 的 叶绿素 相对 含量 ,在 小 麦 09、 玉 米 0、 
棉花 39、 水 稻 59 等 作物 的 氮 素 诊断 与 推荐 施肥 上 得 到 应 用 。 本 研究 将 不 同 航拍 高 度 下 冬小麦 各 样 区 的 色彩 
参数 与 同步 测量 的 SPAD 数值 进行 相关 性 分 析 ( 图 3)。 如 图 3 所 示 ，SPAD 值 与 各 个 备 选 色彩 参数 均 呈 明显 
的 线性 相关 性 ， 除 与 RMR+G+B) 为 负 相 关外 , 与 其 他 3 个 备 选 色彩 参数 均 为 正 相 关 性 , 经 检验 ,相关 性 均 在 
P<0.01 水 平 显著 。SPAD 与 RA(R+G+B)、G/(R+G+B)、B/(R+G+B) 以 及 VARI 的 相关 决定 系数 (R”) 分 别 是 
0.926 9、0.306 2、0.792 9 和 0.898 9。SPAD 与 RA(R+G+B) 相 关 性 最 好 ,其 次 为 VARI。 


0.5 


= 


色彩 参数 Color parameter 


xR/(R+g+b)  。G/(R+G+b) ©¢B/(R+gtb) AVARI 


33 35 37 39 41 43 45 47 49 51 53 
SPAD 


图 3 各 样 区 冬小麦 冠 层 图 像 色彩 参数 与 叶片 SPAD 的 相关 性 


Fig.3 Correlation between canopy photograph color parameters and leaves SPAD of winter wheat 


SPAD 数值 易 受 到 多 种 因素 的 影响 ， 且 存在 饱和 现象 "难以 实现 对 作物 氮 素 营养 的 高 精度 诊断 。 因 此， 
本 研究 采用 冬小麦 茎 基部 硝酸 盐 浓 度 与 各 备 选 色 彩 参 数 做 进一步 相关 性 分 析 ， 以 筛选 出 诊断 冬小麦 氮 素 营 
养 状态 的 最 敏感 的 色彩 参数 。 如 图 4 所 示 , 不 同 航 拍 高 度 获得 的 冬小麦 冠 层 色 彩 参 数 与 冬小麦 茎 基部 硝酸 
盐 含 量 均 呈 窜 函 数 相关 关系 , 旦 R/(R+G+B) 相 关 性 最 好 ,决定 系数 为 0.78， 其 次 为 B/(R+G+B) 和 VARI, 其 
决定 系数 均 为 0.76。 参 考 各 备 选 参 数 与 SPAD 的 相关 性 分 析 ， 可 以 认为 RMR+G+B) 是 反应 作物 氮 素 状态 的 
最 敏感 参数 ， 其 次 为 VARI。R/(R+G+B) 的 取 值 范围 是 从 0.27 到 0.38( 表 1), 数值 变 程 为 0.11， 而 VARI 的 取 
值 范围 从 0.01 到 0.36, 数值 变 程 为 0.35。 从 反 演 作物 氮 素 模型 的 鲁 棒 性 考虑 , R/(R+G+B) 的 变异 系数 为 0.11， 
数据 离散 度 太 低 , 不 确定 性 较 高 ， 容 易 导 致 异常 值 的 出 现 ,而 VARI 的 变异 系数 为 0.39, 不 确定 性 要 远 低 于 
R/(R+G+B)。 在 与 冬小麦 叶片 SPAD 及 茎 基部 硝酸 盐 含量 的 相关 决定 系数 上 ，VARI 略 低 于 R/(R+G+B), 综 
合 考虑 色彩 参数 的 不 确定 性 , VARI 更 适合 作为 监测 冬小麦 氮 素 营养 状态 的 色彩 参数 。 
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图 4 各 样 区 冬小麦 冠 层 图 像 色 彩 参数 与 冬小麦 茎 基部 硝酸 盐 浓 度 的 相关 性 


Fig.4 Correlation between canopy photograph color parameters and stem sap nitrate concentration of winter wheat 


表 1 冬小麦 茎 基部 硝酸 盐 浓度 与 各 色彩 参数 统计 分 析 


Table 1 Statistic analysis of stem sap nitrate concentration and color parameters of winter wheat and their correlations 


色彩 参数 惨 最 小 值 最 大 值 变 幅 
Color parameter Minimum Maximum Amplitude 
R/(R+G+B) 0.80” 0.27 0.38 0.11 
G/(R+G+B) 0.22™ 0.38 0.45 0.07 
B/(R+G+B) 0.74™ 0.21 0.3 0.10 
VARI 0.78™ 0.01 0.36 0.35 


** 表 示 在 P<0.01 水 平 显著 相关 。*#*# means significant relationship at P< 0.01. 


2.2.2 ”适宜 航拍 高 度 的 确定 


无 人 机 在 不 同 高 度 对 冬小麦 进行 航拍 时 受 图 像 分 辨 率 、 观 测 角度 等 因素 的 影响 , 会 导致 图 像 色 
的 变化 。 为 研究 不 同 高 度 航 拍 图 像 对 冬小麦 氮 素 营养 状态 监测 的 差异 , 在 其 他 条 件 (如 相机 设置 、 无 


行 速度 等 ) 均 相同 的 条 件 下 , 分 别 在 8 m、10 m、12 m、14 m、16 m 和 55 m 高 度 对 冬小麦 进行 航拍 ,然后 


提取 各 样 区 的 VARI 数值 与 同步 观测 的 冬小麦 茎 基部 硝酸 盐 浓 度 进行 相关 性 分 析 , 结果 见 图 5, 冬 小 
部 硝酸 盐 浓度 与 不 同 高 度 航拍 图 像 提 取出 的 VARI 均 呈 明显 的 指数 相关 关系 ,并 且 拟 合 的 指数 函数 
为 接近 。 冬 小 麦 茎 基部 硝酸 盐 浓 度 与 8 m、10 m、12 m、14 m、16 m 和 55 m 处 VARI 的 相关 决定 


别 是 0.84、0.84、0.80、0.84、0.82 和 0.78, 在 8~16 m 航拍 高 度 VARI 对 冬小麦 氮 素 营养 的 反应 差异 较 小 ,而 


彩 参 数 
人 机 飞 


麦 茎 基 
曲线 较 
系数 分 


在 55 m 高 度 二 者 的 相关 性 明显 降低 。 因 此, 8~16 m 均 可 作为 冬小麦 航拍 氮 素 诊断 的 适宜 高 度 。 从 航 
的 安全 性 考虑 ， 建议 采 用 16 m 的 航拍 高 度 ,这样 可 以 避 开 地 面 树 木 、 电 杆 等 物体 的 影响 ,方便 作业 
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图 5 不 同 航拍 高 度 航 拍 冬小麦 冠 层 图 像 色彩 参数 VARI 与 茎 基部 硝酸 盐 含量 相关 性 的 比较 


3 作业 


o 


Fig. 5 Comparison of relationships between stem sap nitrate concentration and VARI of aerial photographs of winter wheat 


obtained at different heights 
2.2.3 ”冬小麦 氮 素 诊断 指标 确定 
鉴于 8~16mm 均 适宜 作为 冬小麦 的 航拍 高 度 , 各 航拍 高 度 下 VARI 数 值 相 近 且 与 冬小麦 茎 基部 硝 
度 都 具有 较 好 的 相关 性 ,本 研究 将 8~16 m 航拍 高 度 下 获得 的 各 样 区 的 VARI 数值 进行 平均 ,然后 与 
测 的 冬小麦 茎 基部 硝酸 盐 浓 度 进行 相关 性 回归 建 模 ， 以 获得 利用 航拍 图 像 反 演 冬 小 麦 氮 素 营养 状态 


酸 盐 浓 


同步 观 
的 诊断 


模型 ， 如 图 6 所 示 。VARI 平 均值 与 冬小麦 茎 基部 硝酸 盐 浓度 具有 较 好 的 指数 相关 性 ， 据 此 获得 冬小麦 航拍 


氮 素 营养 诊断 方程 : 冬小麦 茎 基部 硝酸 盐 浓 =2.103 4e™™，。 
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y=2.103 4e18.874> 
R=0.823 1 P<0.01 
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图 6 不 同 高 度 (8~16 m) 航 拍 冬小麦 冠 层 图 像 VARI 与 茎 基部 硝酸 盐 浓 度 的 指数 模型 


Fig.6 Exponential model build by average VARI of aerial photographs from different aerial heights and stem sap concentration of 


结合 


据 VARI 指标 计算 的 不 同 产量 目标 下 冬小麦 的 氮肥 妃 肥 量 。 


标 产 量 时 


winter wheat 


像 的 VARI 反 演 获得 冬小麦 茎 基部 硝酸 盐 浓度 后 ， ee 
量 , 根据 氮 素 平衡 原理 计算 不 同 VARI 对 应 的 氮肥 追肥 量 人 


， 实 现 氮肥 的 精准 管理 。 表 2 为 依 


表 2 


相 | 


不 同 冠 层 图 像 VARI 及 产量 目标 下 冬小麦 氮肥 追肥 量 参 考 指标 


Table 2 Recommended winter wheat topdressing amount of nitrogen fertilizer under different VARIs of aerial photographs 


and yield targets of winter wheat 


VARI <0.2901 0.2901~0.3116 0.3116~0.3269 0.3269~0.3387 0.3387~0.3484 >0.3484 

营养 状况 极 低 低 中 高 较 高 极 高 
Nitrogen nutrition status Very low Low Moderate High Higher Very high 
目标 产量 Yield target (thm2) 推荐 追 氮 量 Recommended topdressing amount of nitrogen fertilizer [kg(N) -hm a 

7.50~8.25 178 139~178 99~139 60~99 21~60 21 
6.75~7.50 162 126~162 90~126 55~90 19~55 19 
6.00~6.75 153 119~153 84~119 50~84 17~50 17 
5.25~6.00 146 114~146 81~114 50~81 15~50 15 
4.50~5.25 104 80~104 56~80 31~56 7~31 7 


2.3” 夏 玉米 氮 素 营养 诊断 分 析 


在 夏 玉米 大 喇叭 


口 期 选择 10 m、15 m、20 m 和 50 m 分 别 进行 航拍 , 研究 不 同 高 度 下 各 备 选 色彩 参数 
对 夏 玉米 氮 素 营养 状态 反映 的 差异 ， 以 确定 敏感 的 色彩 参 


数 及 适宜 的 航拍 高 度 。 同 样 采 用 掩 膜 的 方法 吻 除 


非 叶片 像素 , 统计 出 不 同 航 拍 高 度 下 各 个 样 区 的 R/(R+G+B)、G/(R+G+B)、B/(R+G+B) 以 及 VARI, 与 同步 


测量 的 夏 玉米 叶脉 硝酸 盐 浓 度 进行 相关 性 分 析 ， 以 确定 最 为 敏感 的 色彩 参数 及 适宜 的 航拍 高 度 ,统计 数据 
见 表 3。 
表 3 不 同 航 拍 高 度 下 夏 玉米 冠 层 图 像 色彩 参数 与 叶脉 硝酸 盐 浓 度 的 关系 
Table 1 Relationship analysis between color parameters of aerial photographs obtained at different heights and vein nitrate 
concentration of summer corn 
航拍 高 度 R/(R+G+B) G/(R+G+B) B/(R+G+B) VARI 
Aerial 
photography 拟 合 方程 决定 系数 拟 合 方程 决定 系数 拟 合 方程 决定 系数 拟 合 方程 决定 系数 
height (m) Fitted equation R’ Fitted equation R’ Fitted equation R’ Fitted equation R’ 
10 ]=0.000 Sx 0.14 J=7.4-10 x 0.09 J 王 61.664e? 0.02 y=1.67:10°x 0.30* 
15 =1.28:10*x2137 0.24* y=1.39-10 x 0.02 J]=5.69-10° x 0.07 y=2.02:107x 0.35* 
20 J=1.6:10 1 3 0.42** y=1.37:10 x"! 0.01 y=1.34-100x127 0.18 y=1.49:101x1 0.59** 


刊 


aN i a 
CNINAMIVS 


50 =3.6:10 .x 0.67** y=4.32-10™ x 0.24** y=1.526-10*x3%4 0.92** y=3.90-10°x5% 0.50** 


”+ 是 色彩 参数 ,y 是 夏 玉米 叶脉 硝酸 盐 浓 度 [mg(NOJL ]。* 和 ** 分 别 表示 在 P<0.05 和 P<0.01 水 平 显著 相关 。xis color parameter, y is vein nitrate 
concentration [mg(NO3):L !]. * and ** mean significantly relationship at P < 0.05 and P < 0.01, respectively. 

统计 分 析 显 示 不 同 航拍 高 度 下 同一 个 色彩 参数 与 夏 玉米 叶脉 硝酸 盐 浓 度 的 相关 性 存在 较 大 差异 ,这 意 
味 着 从 不 同 高 度 下 拍摄 的 图 像 其 色彩 参数 存在 变化 ,引起 的 原因 可 能 是 图 像 分 辨 率 以 及 拍摄 角度 不 同 捕获 
的 夏 玉米 冠 层 的 变化 , 具体 原因 需要 做 进一步 的 分 析 , 在 诸多 拟 合 方程 中 , 50 m 高 度 航拍 图 像 的 B/(R+G+B) 
与 夏 玉米 叶脉 硝酸 盐 浓 度 拟 合 方程 的 决定 系数 最 高 (0.92), 在 P<0.01 水 平 下 显著 相关 。 因 此 , 选择 50 m 作 
为 夏 玉 米 的 适宜 航拍 高 度 ,以 B/(R+G+B) 作 为 反映 夏 玉 米 氮 素 营 养 状态 的 最 敏感 参数 , 诊断 方程 为 : 夏 玉 
米 叶脉 硝酸 盐 浓 度 =1.526:10”[B/(R+G+B)J”“*。 利 用 航拍 图 像 反 演 获得 夏 玉米 叶脉 硝酸 盐 浓度 后 ,参考 夏 
玉米 常规 氮肥 追肥 诊断 方法 中 ,可 以 实现 依据 航拍 图 像 B/(R+G+B) 参 数 进行 夏 玉 米 氮 肥 追 肥 量 的 估算 ( 见 
表 4)， 实 现 夏 玉米 氮肥 追肥 的 精准 管理 。 

表 4 不 同 夏 玉 米 冠 层 图 像 色 彩 参数 [B/(R+G+B)] 下 氮肥 追肥 量 参考 指标 (目标 产量 为 7~8 thm-) 
Table 3 Recommended topdressing amount of nitrogen fertilizer under different B/(R+G+B) of aerial photographs obtained at 


50 m high of summer corn with yield target of 7~8 thm2 


0.2606~0.262 0.2626~0.263 0.2639~0.265 0.2650~0.265 0.2658~0.266 0.2666~0.267 
B/(R+G+B) <0.2606 >0.2672 
6 9 0 8 6 2 


推荐 施 氮 量 
Recommended arhount 250 200~250 150~200 100~150 60~100 30~60 0-30 0 
of nitrogen fertilizer 


[kg(N) hm] 


2.4 诊断 指标 验证 

为 验证 诊断 模型 的 精度 , 2016 年 4 月 8 日 和 7 月 22 日 在 肥料 定位 试验 地 对 冬小麦 和 夏 玉米 分 别 以 16 m 
和 50 m 的 高 度 进行 航拍 , 利用 冠 层 图 像 对 作物 氮 素 营养 状态 进行 反 演 , 同时 测量 相同 样 点 冬小麦 茎 基部 和 
夏 玉 米 叶 脉 的 硝酸 盐 浓 度 。 冬小麦 和 夏 玉 米 氮 素 营养 状态 的 实测 值 与 预测 值 比较 如 图 7 所 示 。2016 年 利用 
航拍 图 像 反 演 的 冬小麦 和 夏 玉 米 氮 素 状 态 与 实测 数据 具有 很 好 的 拟 合 精度 ， 反 演 与 实测 的 冬小麦 葵 基 部 硝 
酸 盐 浓度 决定 系数 为 0.80, 在 P<0.01 水 平 显 赣 相 关 ; 夏 玉 米 叶 肪 硝酸 盐 浓度 的 诀 定 系数 为 0.85， 同 样 在 
P<0.01 水 平 显著 相关 。 结 果 表 明 本 研究 所 建立 的 冬小麦 和 夏 玉 米 氮 素 营养 诊断 指标 可 行 。 
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图 7 冬小麦 、 夏 玉米 氮 素 营养 航拍 诊断 精度 


Fig.7 Verification of nitrogen nutrition diagnosis using aerial photographs of winter wheat and summer corn 


2.5 无 人 机 航拍 氮 素 营养 诊断 应 用 
利用 上 述 研究 建立 的 航拍 方法 及 诊断 模型 ，2017 年 在 河北 省 沧州 市 南 皮 县 东 五 拨 村 (38.068*N， 
116.825°E) 约 30 hm 的 平整 农田 进行 了 航拍 氮 素 营养 诊断 应 用 。 选 择 在 冬小麦 拔节 期 和 夏 玉 米 大 喇叭 口 期 
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进行 航拍 作业 , 通过 软件 开发 建立 了 图 像 诊 断 批 处 理 系 统 ， 实现 航拍 图 像 经 纬度 、 叶 片 掩 膜 、 色彩 参数 以 
及 氮肥 追肥 量 的 快速 批量 处 理 , 满足 农田 管理 需求 。 利 用 航拍 诊断 的 结果 ,通过 反 距离 加 权 插 值 法 获得 每 
个 田 块 的 氮肥 追肥 量 ,为 便于 作业 ， 以 尿素 为 例 , 给 出 了 氮肥 追肥 作业 图 。 图 8 中 白 点 为 航拍 图 像 的 位 置 轨 
迹 , 长 条 为 地 块 。 经 过 地 面 调查 和 农户 走访 , 作业 图 给 出 的 氮肥 追肥 量 高 低 与 作物 长 势 强 弱 在 空间 分 布 上 
较为 吻合 。 
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图 8 基于 航拍 图 像 诊 断 的 冬小麦 (a)、 夏 玉米 (b) 氮 肥 追 肥 施 肥 量 [kg( 尿 素 )-hm] 
Fig. 8 Prescription maps of top-dressing nitrogenous fertilizer amounts [kgCurea)hm2] for winter wheat (a) and summer corn (b) by 


aerial photograph diagnosis 


3 讨论 与 结论 
解决 农业 生产 对 化 肥 的 依赖 以 及 由 此 造成 的 面 源 污染 等 问题 , 实现 化 肥 施 用 量 零增长 的 目标 , 需要 切 
实 可 行 的 作物 养分 诊断 与 施肥 技术 , 特别 是 对 氮 素 的 快速 、 无 损 、 低 成 本 诊断 技术 。 叶 色 诊 断 方 法 简单 易 
了 , 并且 一 直 以 来 就 被 我 国 农民 所 采用 , 但 过 于 定性 且 依 靠 个 人 经 验 。 数 码 相 机 的 采用 使 得 可 以 应 用 数字 
图 像 技 术 , 定量 分 析 作 物 叶 片 对 红 光 、 蓝 光 和 绿 光 的 吸收 差异 , 进而 追 湖 引起 该 差异 的 氮 素 丰 缺 原因 ， 实 现 
氮 素 营养 无 损 、 低 成 本 的 定量 化 分 析 和 决策 户 。 随 着 士 地 流转 和 农业 适度 规模 经 营 已 成 为 趋势 ， 必 须 满足 
大 面积 生产 对 氮 素 诊断 的 需求 ,即使 采用 数字 图 像 诊 断 方法 ,拍照 的 工作 量 也 是 很 大 , 特别 是 玉米 等 高 秆 
作物 ， 难 以 进入 农田 进行 大 量 拍照 。 

面 对 新 的 生产 需求 ， 本 研究 基于 数字 图 像 诊断 方法 的 研究 基础 ， 利 用 无 人 机 航拍 技术 ,自制 了 数码 相 
机 搭载 平台 , 通过 航线 制定 和 自动 航拍 , 解决 了 作物 冠 层 图 像 快速 获取 问题 。 针对 冬小麦 和 夏 玉 米 , 在 适宜 
航拍 高 度 以 及 对 作物 氮 素 营养 敏感 色彩 参数 的 筛选 上 进行 了 深入 研究 。 结 果 表 明 : 在 8~16 m 高 度 冬小麦 航 
拍 图 像 的 VARI 对 其 氮 素 营养 均 很 敏感 ， 考虑 到 航拍 安全 ,建议 冬小麦 采用 16 m 的 航拍 高 度 ; 不 同 航拍 高 
度 下 夏 玉 米 图 像 的 色彩 参数 对 其 氮 素 营养 的 敏感 性 存在 较 大 差异 , 50 m 高 度 航拍 图 像 的 B/(R+G+B) 与 夏 玉 
米 叶 脉 硝 酸 盐 浓度 相关 性 最 好 ， 可 以 作为 航拍 诊断 的 依据 。 基于 上 述 分 析 , 2016 年 进行 了 验证 , 结果 表明 航 
拍 诊断 模型 在 冬小麦 和 夏 玉 米 上 均 有 较 好 的 监测 精度 。 

相对 于 传统 氮 素 化 验 分 析 方 法 ,数字 图 像 诊断 法 更 为 便捷 ,但 要 实现 在 生产 中 的 成 熟 应 用 ,还 要 做 
一 步 的 深入 研究 和 技术 完善 。 就 本 研究 而 言 , 仅 对 长 期 肥料 定位 试验 下 的 单一 品种 进行 了 分 析 , 品种 间 的 
叶 色 差异 需要 考虑 ; 本 文采 用 比值 型 色彩 参数 来 消除 光线 对 成 像 色 彩 的 影响 ， 王 娟 等 已 利用 灰 板 校正 以 消 
除 棉花 不 同 生 育 期 图 像 颜色 特征 值 的 亮度 差异 来 提高 诊断 植株 含水 量 的 精度 ， 该 方法 可 以 参考 。 另 外 , 干 
旱 能 够 引起 作物 时 色 和 植株 硝酸 盐 浓度 的 变化 , 这 将 如 何 影响 航拍 图 像 诊断 精度 ,色彩 参数 是 否 存在 饱和 
现象 等 , 这 些 问 题 将 作为 本 研究 的 后 续 工作 的 重点 内 容 ,以 进一步 完善 无 人 机 航拍 氮 素 诊断 方法 ， 并 在 其 
他 作物 上 得 到 应 用 , 促进 精准 施肥 工作 的 开展 。 
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